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To  capture  global  responses  to  metal  poisoning  and  mechanistic  insights  into  metal  toxicity,  gene  expression 
changes  were  evaluated  in  whole  adult  male  zebrafish  following  acute  24  h  high  dose  exposure  to  three  metals 
with  known  human  health  risks.  Male  adult  zebrafish  were  exposed  to  nickel  chloride,  cobalt  chloride  or  sodium 
dichromate  at  concentrations  corresponding  to  their  respective  96  h  LC20,  LC40  and  LCeo  (i.e.  96  h  concentrations 
at  which  20%,  40%  and  60%  lethality  is  expected,  respectively).  Histopathology  was  performed  on  a  subset  of 
metal-exposed  zebrafish  to  phenotypically  anchor  transcriptional  changes  associated  with  each  metal  exposure. 
Here  we  describe  in  detail  the  contents  and  quality  controls  for  the  gene  expression  and  other  data  associated 
with  the  study  published  by  Hussainzada  and  colleagues  in  BMC  Pharmacology  and  Toxicology  (Hussainzada 
et  al.,  2014)  with  the  data  uploaded  to  Gene  Expression  Omnibus  (accession  number  GSE50648). 

©  2014  The  Authors.  Published  by  Elsevier  Inc.  This  is  an  open  access  article  under  the  CC  BY-NC-ND  license 

(http://creativecommons.0rg/licenses/by-nc-nd/4.O/). 


Specifications 

Organism/celi  iine/tissue 

Danio  rerio  (Tubingen  strain) 

Sex 

Male 

Sequencer  or  array  type 

Agilent  4  x  44  K  custom  design  (see  Microarray 
hybridization  for  details) 

Data  format 

Raw  data:  .gpr  files;  normalized  data:  SOFT,  MINML, 
and  TXT 

Experimental  factors 

Metal  exposure  vs.  unexposed  normal 

Experimental  features 

Male  adult  zebrafish  were  exposed  to  nickel  chloride, 
cobalt  chloride  or  sodium  dichromate  concentrations 
corresponding  to  their  respective  96  h  LC20,  LC40  and 
LCeo-  Histopathology  was  performed  on  a  subset  of 
metal-exposed  zebrafish  to  phenotypically  anchor 
transcriptional  changes  associated  with  each  metal  as 
determined  by  microarray  analysis. 

Consent 

None 

Sample  source  location 

In-house  aquaculture  at  US  Army  Center  for 
Environmental  Health  Research,  Eort  Detrick, 

Erederick,  MD,  USA 

*  Corresponding  author  at:  Excet,  Inc.  USACEHR  568  Doughten  Drive  Et.  Detrick,  MD 
21702-5010,  USA. 


Direct  link  to  deposited  data 

Deposited  data  can  be  found  here:  http://www.ncbi.nlm.nih.gov/ 
geo/query/acc.cgi?acc=GSE50648. 

Experimental  design,  materials  and  methods 

Experimental  overview 

Fig.  1  shows  the  design  scheme. 

Fish  exposures 

Briefly,  adult,  male  zebrafish  were  exposed  to  concentrations  of  each 
metal  necessary  for  20%  (LC20;  low),  40%  (LC40;  mid)  and  60%  (LCeo; 
high)  lethality  based  on  prior  96  h  range-finding  experiments.  Expo¬ 
sures  were  conducted  for  24  h  using  control  (no  toxicant)  plus  the 
high,  mid-  and  low  concentrations  of  each  metal  (Table  1 ).  The  intent 
was  to  evaluate  levels  of  toxicant  sufficient  to  induce  a  measurable  re¬ 
sponse  without  producing  lethality  at  24  h.  All  metal  concentrations  in 
the  test  tanks  were  verified  both  before  and  after  exposures.  Initially 
25  adult  (6-9  months)  presumptive  male  zebrafish  were  selected  per 
condition  to  ensure  that  20  male  fish  were  available  for  subsequent  mi¬ 
croarray  analysis  and  histopathology.  At  termination  of  each  study,  five 
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Fig.  1. 


fish  per  condition  were  sacrificed  by  terminal  dose  of  MS222  and  then 
immediately  preserved  in  a  modified  Davidson's  solution  for  histologi¬ 
cal  examination.  For  transcriptional  analysis,  the  remaining  20  zebrafish 
were  sacrificed  by  terminal  dose  of  MS222  and  then  immediately  im¬ 
mersed  whole  in  liquid  nitrogen  and  stored  at  —80  °C  until  RNA 
processing. 


RNA  processing 

Flash  frozen  whole  fish  were  individually  pulverized  under  liquid  ni¬ 
trogen  using  a  SPEX  6750  freezer  mill  (SPEX  Sample  Prep,  Metuchen, 
NJ).  Six  milliliters  of  Trizol®  (Invitrogen,  Carlsbad,  CA)  was  added  to 
the  pulverized  material,  vortexed  well  and  transferred  to  a  10  ml 
Dounce  homogenizer.  An  additional  2  ml  of  Trizol®  was  added  to  re¬ 
move  any  remaining  pulverized  material  from  the  initial  tube  and  trans¬ 
ferred  to  the  same  Dounce.  The  pulverized  material  was  further 
homogenized  by  hand  with  10  strokes  using  Pestle  A  followed  by  10 
strokes  using  Pestle  B.  After  homogenization,  the  manufacturer's  sug¬ 
gested  protocol  was  followed  for  RNA  extraction  with  an  extra  clarifica¬ 
tion  centrifugation  step  to  remove  bone,  scales,  and  other  insoluble 
debris  from  the  homogenate.  Total  RNA  isolated  via  Trizol®  was  then 
subjected  to  further  purification  using  RNeasy®  Midi  kits  (Qiagen, 
GmbH,  Germany)  to  remove  any  residual  salts  or  organic  solvents.  Qual¬ 
ity  and  quantity  of  isolated  total  RNA  were  analyzed  using  an  Agilent 
Bioanalyzer  2100  (Agilent,  Santa  Clara,  CA)  and  verified  with  the 
NanoDrop  ND-1000  Spectrophotometer  (NanoDrop,  Wilmington,  DE). 


Table  1 

Nominal  and  actual  concentrations  (mg/L)  of  metals  used  in  24- hour  zebrafish  exposures. 
From  Hussainzada  et  al.  (2014). 


Treatment 

Nominal 

LC20 

(low) 

Measured 

LC20 

(low) 

Nominal 

LC40 

(mid) 

Measured 

LC40 

(mid) 

Nominal 

LCeo 

(high) 

Measured 

LCeo 

(high) 

NiCb 

45 

42.4 

54 

51.0 

62 

64.0 

C0CI3 

39 

39.7 

50 

46.3 

65 

59.5 

Na2Cr207 

53 

56.5 

65 

69.9 

76 

80.6 

Total  RNA  showing  signs  of  degradation  or  a  RIN  value  <8  was  not  in¬ 
cluded  in  further  analysis. 

Screening 

Prior  to  performing  microarray  analysis,  the  RNA  was  screened 
against  a  primer  panel  to  verify  that  the  RNA  was  isolated  from  male 
fish  only.  An  aliquot  of  each  total  RNA  sample  was  transcribed  to 
cDNA  using  the  Advantage®  RT-for-PCR  Kit  (Clontech,  Mountain 
View,  CA).  Amplification  and  detection  were  performed  on  an  MJ 
Opticon®  2  (BioRad,  Hercules,  CA)  using  SYBR®  Green  PCR  Master 
Mix  (Applied  Biosystems,  Carlsbad,  CA)  with  primers  designed  to  detect 
transcripts  coding  for  vitellogenin  1  {vitl,  expressed  only  in  female  liver) 
and  glyceraldehyde  3-phosphate  dehydrogenase  (gapdh).  gapdh  was 
used  as  an  internal  control  for  normalizing  sample  concentrations. 
Any  sample  which  showed  the  presence  of  female  RNA  (i.e.  high  levels 
of  vitl )  was  eliminated  from  further  analysis. 

Microarray  hybridization 

In  order  to  maximize  statistical  power  while  minimizing  cost,  equal 
amounts  of  total  RNA  from  four  or  five  individual  fish  from  each  expo¬ 
sure  condition  in  each  of  the  four  exposure  experiments  were  combined 
to  create  four  biological  replicate  pools.  Each  biological  replicate  pool 
was  hybridized  to  a  separate  microarray  for  a  total  of  16  microarrays 
(i.e.  four  control  replicates,  four  low  dose  replicates,  four  mid-dose  rep¬ 
licates,  and  four  high  dose  replicates). 

Eor  this  study,  custom  designed  microarrays  were  used.  The  arrays 
were  designed  in-house  using  the  eArray  microarray  design  tool 
(https://earray.chem.agilent.com/earray/;  Agilent  Technologies,  Inc.) 
using  zebrafish  gene  targets  derived  from  Ensembl  build  46  (Zv7  ge¬ 
nome  build)  and  Vega  build  26.  Each  array  contains  44,000  60-mer  oli¬ 
gonucleotides  representing  21,904  zebrafish  gene  targets  with  two 
probes  per  transcript  wherever  possible;  only  94  transcripts  had  a  single 
probe.  Probes  were  designed  using  only  genes  which  are  annotated  in 
published  databases.  All  arrays  were  manufactured  by  Agilent. 

All  microarrays  were  processed  following  Agilent's  One-Color 
Microarray-Based  Gene  Expression  Analysis  Protocol  (Version  5.5, 
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February  2007)  for  processing  4  x  44  K  array  slides.  An  initial  1  |ug  pooled 
total  RNA  input  was  used  with  an  18  hour  overnight  hybridization  at 
65  °C  A  final  wash  using  the  Stabilization  and  Drying  solution  (Agilent 
Technologies,  Inc.)  was  included  after  the  required  specificity  washes  in 
order  to  prevent  any  ozone  related  degradation  of  signal  prior  to  scanning 
the  arrays.  Arrays  were  scanned  on  a  GenePix  Autoloader  4200  AL 
(Molecular  Devices,  Union  City,  CA)  using  the  532  nm  (green)  laser  at 
45%  with  the  green  filter,  PMT  =  400,  and  the  scan  resolution  =  5  jum. 
Raw  images  were  processed  using  GenePix  Pro  6.0  (Molecular  Devices) 
with  the  following  settings:  GAL  file  =  0177662_D_20070905_ 
updateS.gal;  local  background  subtraction;  background  subtraction: 
width  of  background  =  2  feature  diameters;  resize  features:  minimum 
diameter  =  80%;  resize  features:  maximum  diameter  =  110%;  feature 
movement:  maximum  translation  =  20  jum;  and  CPI  threshold  =  0. 
These  setting  were  determined  to  be  optimal  for  grid  alignment  based 
on  optimization  experiments  performed  in  house.  Prior  to  performing 
final  quality  control  checks,  arrays  were  visually  inspected  to  verify  grid 
alignment  and  the  absence  of  any  defects.  All  microarray  data  from  this 
study  have  been  deposited  in  NCBFs  Gene  Expression  Omnibus  [1] 
under  the  accession  number  GSE50648. 

Quality  control  &  normalization 

GenePix  Pro  offers  intensity  data  in  several  alternate  forms.  Eor  this 
study  [2],  median  signal  intensity  and  signal-to-noise  ratio  (SNR)  were 
used.  Median  signal  intensity  is  calculated  as  the  median  value  of  all 
the  pixels  that  are  completely  with  the  feature  indicator;  any  pixel 
touching  the  boundary  is  excluded.  This  software  calculates  SNR  as 
the  difference  between  median  spot  intensity  signal  and  median  back¬ 
ground  divided  by  the  standard  deviation  of  the  background  signal. 

Median  signal  intensity  array  data  was  imported  into  Partek 
Genomics  Suite  (version  6.6;  Partek  Inc.,  St.  Louis,  MO)  for  additional 
analysis.  To  verify  performance  of  the  replicates,  Pearson's  R^  was  calcu¬ 
lated  for  each  set  of  replicates  within  each  toxicant  group.  Any  arrays 
with  R^  <  0.95  were  excluded  from  analysis.  All  arrays  included  in  the 
final  analysis  had  R^  >  0.97  with  their  replicates. 

To  further  increase  confidence,  only  unsaturated  probes  (i.e.  probes 
with  signal  intensity  less  than  65,535)  with  a  SNR  greater  than  or  equal 
to  three  (SNR  >  3)  were  selected  for  analysis.  The  SNR  data  was 
imported  into  Partek  Genomics  Suite;  a  transposed  spreadsheet  was 
created  and  then  transferred  to  Microsoft  Excel  (Microsoft  Corporation, 
Redmond,  WA)  where  a  list  was  created  of  all  probe  sets  with  a  SNR  >  3 
in  all  replicates  of  at  least  one  condition.  This  list  was  then  used  to  filter 
the  median  signal  intensity  data  in  Partek.  Only  unsaturated  probes 
with  SNR  >  3  were  selected  for  analysis  which  yielded  a  subset  of 
15,818  probes  which  map  to  7909  genes.  Quantile  normalization  across 
arrays  was  performed  to  control  for  any  inter-array  variability  then 
these  normalized  probe  intensities  were  log  transformed. 

Basic  analysis 

Using  Partek  Genomics  Suite,  three  sets  of  ANOVAs  (one  set  for  each 
metal)  were  performed  to  identify  probes  that  were  differentially 
expressed  between  each  metal  exposure  group  and  its  control  group. 
The  ANOVAs  included  terms  for  treatment  (exposed  or  unexposed) 
and  concentration  (control,  low,  mid,  or  high)  and  an  interaction  term 
for  “treatment  *  concentration”.  Eor  each  chemical,  contrasts  were  per¬ 
formed  to  determine  significance  between  each  dose  group  and  its  re¬ 
spective  control.  A  step-up  Benjamin!  and  Hochberg  false  discovery 
rate  (EDR)  of  0.01  was  used  to  select  differentially  expressed  probes. 
An  EDR  a  =  0.01  was  selected  as  a  cut-off  for  the  combined  nickel,  chro¬ 
mium,  and  cadmium  datasets  including  all  replicate  pools.  Probes  not 
meeting  the  EDR  a  =  0.01  threshold  were  eliminated  from  further  anal¬ 
ysis.  The  resulting  list  was  further  refined  by  submitting  it  to  a  second 
filter  which  specified  a  1.8-fold  change  between  exposed  versus  control. 
Only  transcripts  with  probes  which  passed  both  filters  were  included  in 


Fig.  2.  Total  number  of  genes  significantly  (FDR  =  0.01)  altered  by  at  least  1.8  fold  in 
response  to  chromium  (Cr),  cobalt  (Co),  or  nickel  (Ni)  versus  untreated  controls. 

From  Hussainzada  et  al.  (2014). 

the  final  analysis.  Eor  all  transcripts  with  two  probes,  fold  changes  for 
each  probe  pair  were  averaged  to  generate  a  single  value  for  each  tran¬ 
script.  Eor  all  transcripts  with  a  single  probe,  fold  change  data  stood  as  is. 

Exposure  to  chromium,  cobalt,  or  nickel  significantly  altered  the 
expression  of  696, 461,  and  287  genes  respectively  (Eig.  2). 

Discussion 

We  describe  here  a  unique  dataset  of  male  zebrafish  transcriptomic 
responses  to  acute  metal  toxicity.  This  dataset  is  composed  of  genome¬ 
wide  gene  expression  data  measured  using  a  custom  Agilent  array  plat¬ 
form.  The  gene  expression  data  is  of  high  quality  and  consistent  with 
histopathologic  data  collected  during  this  study  [2].  Using  this  data, 
we  identify  changes  in  gene  expression  of  groups  of  genes  consistent 
with  adaptive  responses  induced  by  exposure  to  nickel,  cobalt,  or  chro¬ 
mium  including  acute  phase  response,  cell  cycle  regulation,  apoptosis, 
and  metabolic  depression. 

Disclaimer 

Opinions,  interpretations,  conclusions,  and  recommendations  are 
those  of  the  authors  and  are  not  necessarily  endorsed  by  the  U.S. 
Army.  Research  was  conducted  in  compliance  with  the  Animal  Welfare 
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Guide  for  the  Care  and  Use  of  Laboratory  Animals  (NRC  201 1 )  in  facili¬ 
ties  that  are  fully  accredited  by  the  Association  for  the  Assessment  and 
Accreditation  of  Laboratory  Animal  Care,  International.  The  research  de¬ 
scribed  herein  was  sponsored  by  the  U.S.  Army  Medical  Research  and 
Materiel  Command,  Military  Operational  Medicine  Research  Program. 
Citations  of  commercial  organizations  or  trade  names  in  this  report  do 
not  constitute  an  official  Department  of  the  Army  endorsement  or 
approval  of  the  products  or  services  of  these  organizations. 

This  project  was  supported  in  part  by  an  appointment  to  the  Research 
Participation  Program  for  the  U.S.  Army  Medical  Research  and  Materiel 
Command  administered  by  the  Oak  Ridge  Institute  for  Science  and 
Education  through  an  agreement  between  the  U.S.  Department  of  Energy 
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